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I. INTRODUCTION 



Since Loss and DiVincenzo proposed the utilization of quantum dots as the physical 
implementation of the qubit^, there has been a huge amount of work devoted to tackle 
the numerous and subtle difficulties involved in the problem. There are some excellent 
reviews^ and books^ that summarize the progress experimented by the field, but it is 
extremely difficult to keep up with the new developments. 

As much as any other proposal to implement a qubit, the spin degree of freedom of an 
electron trapped in a quantum dot (QD), the original proposal made by Loss and DiVincenzo, 
must face a number of challenges owed to the intrinsic physics that governs its behavior as 
a qubit. Not to mention the challenges offered by other physical implementations that try 
to catch the attention of the community^. The double quantum dot scheme^ was devised 
to circumvent the unavoidable decoherence induced by the interaction between the angular 
momentum of the electron with the nuclear spins of the atoms that form the QD^. For this 
scheme, the realization of multiple qubit quantum gates has been showiA Nevertheless, since 
the coupling between single QD's seems a bit problematic and maybe even more involved 
in the double QD scheme, there has been a number of proposal showing that it is possible 
to implement quantum control and refocusing techniques in single quantum dots to restore 
the role of the single QD as a bona fide qubit. This, together with techniques designed to 
distinguish between spatial states of the trapped electron, make interesting again the search 
of new one-electron structures that can be controlled with the exquisite precision required 
for the quantum information tasks. 

The control of quantum systems, at least when the decoherence mechanisms are absent 
or "turned-off " , is implemented using pulses of external fields or manipulating an adequate 
parameter. A complete knowledge of the spectrum allows the use of "navigating methods" 
that make possible going from (almost) any initial state to the desired target stated Never- 
theless, the most used approach to achieve the switching between the two qubit basis states 
is the optimal control theory*^. The application of Krotov's algorithm^ usually leads to 
speedups of the transition time. This method has successfully been applied to one— ^ anc j 
two-electron quantum dots*^, allowing fast charge transfer with larger fidelities than the 
obtained with simpler sinusoidal pulses. As much promising as the optimized pulse method 
seems, the introduction of a complex modulation of the control field necessarily introduces 

2 



a host of new error sources that have not been properly analyzed. In this sense, achieving 
a good control of the switching between states of a quantum dot without resorting to a 
complicated pulse sequence is worth of study, and is the aim of this work. 

Recently, it has been shown that multi-layered quantum dots can be designed to se- 
lectively modulate the spatial extent of the electronic density of its eigenstates^ 1 ^. This 
feature, together with the dipole selection rule permits, as we will show, the design of a 
quantum dot with two states that can be switched, robust and efficiently, using only sinu- 
soidal pulses. These two states are the basis states of our model qubit. Using high-precision 
ab-initio numerical calculations and exact solutions, where available, we aim to study the 
spectrum, eigenstates of the quantum dot and dynamical properties of the qubit. 

The paper is organized as follows, in Section [III a realistic multi-layered quantum dot 
model is presented and qualitatively analyzed. In Section IIHI the properties of the eigen- 
values and eigenstates of the model are obtained. The study of the spatial extent of the 
eigenstates allows the identification of potential qubit basis states. The time evolution of the 
quantum state when the system is driven by an external sinusoidal radio frequency field (rf) 
is the subject of Section IIVI It is shown that the sinusoidal driving is enough to obtain an 
excellent switching between the qubit basis states with very low probability leakage, making 
unnecessary the utilization of complicated envelope functions as is customary in optimal 
control theory. The stability of the qubit oscillations is also tested considering the effect 
of off-resonance driving. In Section [V] we briefly present a different model potential that 
also define an excellent qubit, whose properties are analyzed through the lines drawn in 
Sections [III IHH and [IVl Finally, we discuss our results in Section ED 

II. THE MODEL 

In the Effective Mass Approximation, the Hamiltonian of trapped particles assumes a 
simple form since the many-body interactions are reduced to a bounding potential. In 
a spherical layered quantum dot, where each layer is made of a different material, the 
bounding potential is given by the conduction band off-sets of each material. Figure [T] 
shows, schematically, the radial profile of the bounding potential for a quantum dot made of 
two different materials. The basic structure consists of a central core and two wells separated 
by a barrier—. 
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FIG. 1. (Color on-line) a) The electronic density for a quantum state well localized in the innermost 
potential well of a layered quantum dot. The radial step-like potential given by Equation [1] is also 
shown (red dashed line), b) The electronic density jumps from the innermost potential well to 
the outermost one when the radius of the quantum dot central core is changed. The target-like 
pattern to the right of the figure corresponds to the cross-section of the quantum dot, the grey 
zones correspond to both potential wells, while the central core, and the barrier are depicted with 
white. 



Despite its simplicity, the piecewise potential takes into account a number of experimental 
features and allows the formulation of accurate and simple models for many nano-structures. 
The current semiconductor technology permit the fabrication of layered structures where the 
radius of each layer can be tuned with great precision. Figure [1] also shows the qualitative 
behavior of a given eigenstate, when the central core is wide enough the electronic density is 
mostly located in the inner potential well, conversely, when the radius of the central core is 
diminished the electronic density jumps to the outermost potential well for a certain critical 
value. So, changing the quantum dot architecture is equivalent to choose between different 
spectra and sets of eigenstates, whose physical properties can be dramatically changed just 
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altering the design of the quantum dot. In particular, the change in the spatial extent 
of low lying eigenstates and the modulation of the oscillator strength associated to these 
eigenstates was analyzed in the work by Ferron et al.—. 

More precisely, in this work the bounding potential considered is given by 
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where Vq = 0.9 eV which corresponds to the band offset between CdSe and ZnS, while the 



effective masses are Tn* CdSe = 0.1 3m e and m* ZnS 



0.28m e , m e is the mass of the bare 



electron^. 

Since we are interested in the properties of a one electron quantum dot when it is driven 
by an external field, the Hamiltonian takes the form 



H = H + V ext (r,t) 



(2) 



where 



ft 2 

H n = --V 



V + VJr) 



(3) 



2 \m(r) 

where m(r) is the position-dependent effective mass of the electron, and the bounding poten- 
tial is given in dJ The kinetic energy term in Equation E] preserves the Hermitian character 
of the Hamiltonian operator when the mass is position dependent^. In the present case the 
mass a step-like function of the radial coordinate. 



III. SPECTRUM AND EIGENSTATES 

The spectrum and eigenstates of Hamiltonian [3] can be obtained exactly or numerically. 
In any case, the problem has spherical symmetry so the eigenvalues depend on two quantum 
numbers, (n r ,£), where n r is radial quantum number, or the number of nodes of the radial 
eigenf unction, and £ is the orbital angular momentum quantum number. On the other hand, 
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note that since the exact solution of the eigenvalue problem involves the roots of transcen- 
dental equations, it can be difficult to figure out how many quasi-degenerate eigenvalues has 
the problem. This difficulty is particularly cumbersome for large values of r c . Conversely, 
the variational methods detect fairly well almost degenerate eigenvalues, so the data shown 
in this Section was double-checked comparing the results from the numerical and analytical 
procedures, 

The numerical solutions were obtained using B-splines basis sets, which are well suited to 
implement the boundary conditions imposed by the step-like nature of the potential and the 
effective mass. Besides, B-splines results are very accurate in comparison with calculations 
based on Gaussian, Hylleraas and finite-element basis sets^. 

To use the B-splines basis, the normalized one-electron orbitals are given by 

Mr) = C n ^-^; n = l,... (4) 

where B^+i{r) is a B-splines polynomial of order k. The numerical results are obtained by 
defining a cutoff radius R, and then the interval [0, R] is divided into / equal subintervals. 
B-spline polynomials^ (for a review of applications of B-splines polynomials in atomic and 
molecular physics, see ref.— ) are piecewise polynomials defined by a sequence of knots t± = 
< t<i < ■ ■ ■ < t2k+i-i = R and the recurrence relations 

flf(r) = ( 1 " r < tM . (5) 
I otherwise, 

B« (r) = r ~\ Bt 1] (r) + (r) (* > 1) . (6) 

In this work, we use the standard choice for the knots in atomic physics^ t\ — ■ ■ ■ — tk — 
and tk+i = ■ ■ ■ = t 2 k+i-i = R- We choose an equidistant distribution of inside knots. The 
constant C n in Eq.(3) is a normalization constant obtained from the condition (n\n) = 1, 



1 - ll/2 - 



Because -Bi(O) ^ and Bi + k-i(R) ^ 0, we have N = I + k — 3 orbitals corresponding to 
B 2 , . . . , B I+k _ 2 . In all the calculations we used the value k — 5, and, we do not write the 
index k in the eigenvalues and coefficients. 



6 



I 



2 4 6 8 10 

r c (nm) 



FIG. 2. (Color on-line) The spectrum of a typical layered quantum dot as a function of the inner 
core radius, r c . The radius of each layer is, r\ = r c + 0.8, r 2 = r\ + 3.5, and r 3 = r 2 + 1, all in 
nanometers. The black solid lines correspond (from bottom to top) to eigenvalues with quantum 
numbers (n r = 0,£), £ = 0,1,2,... ,22, while the red solid lines correspond to (from bottom to 
top) eigenvalues with quantum numbers (n r = 1,£), £ = 0,1,2, ... , 12. 

Figure |2] shows the spectrum of a double- well quantum dot as a function of the central 
core radius, the width of the two wells and the barrier are kept constant, so r\ = r c + 0.8nm, 
r 2 — ri+3.5nm and r 3 = r 2 + lnm. As can be appreciated, the spectrum is quite complicated 
and there is not a couple of energy values well separated from the others, which is a common 
criterion to identify possible basis eigenstates for a qubit. As we will show, the availability 
of a couple of eigenstates that can be easily switched with a simple pulse depends not only 
on the characteristics of the spectrum, but also in the eigenstates spatial extent. 

Figure [3] shows the electronic density corresponding to all the bounded eigenstates of two 
similar devices. Panel a) corresponds to a device with only two eigenstates well localized in 
the innermost potential well , while panel b) shows a device with three eigenstates localized 
in the innermost potential well. The two devices only differ in the central core radius size, 
it is bigger for the b) case. The radii of both devices are, Device l:r c = 1 nm, r\ = 1.8 nm, 
r 2 = 5.3 nm and r 3 = 6.3 nm. Device 2: r c = 2 nm, r± = 2.8 nm, r 2 = 6.3 nm and r 3 = 7.3 
nm. 

In both cases the eigenstates localized in the inner potential well have different angular 




FIG. 3. (Color on-line) Electronic densities for the ground and excited states for two different 
CdSe/ZnS devices, (a) Device 1, (b) Device 2. The electronic densities of the eigenstates corre- 
sponding to the qubit states are denoted with q\ and q%, which have quantum numbers (1,0) and 
(1,1), respectively. . The black dashed vertical lines show the positions of the barriers and wells 
of the quantum dot. Both Devices are particular cases of the layered quantum dot whose spectra 
are depicted in Figure [2j 

momentum quantum numbers, in a) A£ = 1. Actually, in the Device 1 case, the qubit 
lower state has quantum numbers n r — 1 and i = 0, while the higher one has n r — 1 and 
£ — 1, while in the Device 2 case the third state has n r — 1 and 1 = 1. Figure [3] shows 
two remarkable facts, i) the number of eigenstates localized in a given potential well can 
be chosen for realistic quantum dots parameters and, ii) the number of states of election is 
pretty stable against small changes in the device's parameters (see also Figure However, 
as we will show latter a change in the number of states localized in the potential well of 
interest can produce a huge change in the dynamical behavior once a external driving is 
applied to the device. 

On the other hand, despite that we mostly present results for particular sets of param- 
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eters (effective masses, radii, etc.) the physical traits that, at some extent, guarantee the 
existence of a small number of eigenstates well localized in a multi-well potential are fairly 
general. Indeed, later on we will show a example with a continuous bounding potential 
whose spectrum, eigenstates and dynamical behavior are strikingly similar to those of the 
devices with step-like potentials. 

IV. SINUSOIDAL DRIVING OF THE ELECTRON 

The most simple non-trivial driving, both from an numerical and experimental point of 
view, that can be applied to the trapped electron is 



where Aq and u are the strength and the frequency of the driving, respectively. The time- 
dependent potential, Equation El models the effect of a time-periodic spatially constant 
electric field applied in the z direction. As the potential Equation depends on z, the 
dipole selection rules impose transitions between eigenstates that differ in orbital angular 
momentum, A£ = ±1. 

The time evolution of the electron quantum state is governed by the Schrodinger equation, 



where \& is the quantum state and H is the Hamiltonian in Equation [2j Since we want to 
analyze the behavior of the devices presented in the previous Sections as qubits, we will 
study the electron quantum state evolution, taking as initial condition the lowest eigenstate 
that is localized in the innermost potential well, from now on the |0) state of our putative 
qubit, this state has i = 0. The other qubit basis state, |1) is the eigenstate with i = 1 that 
is also localized in the innermost potential well. Ideally, to qualify as a qubit, a physical 
system would perfectly switch between the two basis states under the appropriate driving, 
i.e. if \c q i\ 2 and \c q 2\ 2 are the time dependent probabilities that the electron is in the |0) or 
in the |1) state, respectively, then \c ql \ 2 + c q2 \ 2 = 1. 

Except for ideal two level systems, there is a finite probability that after a switching 
operation |c 9 i| 2 + c q2 \ 2 < 1, i-e. the qubit leaks probability. The leakage, defined as 1 — 
(|c g i| 2 + c q2 \ 2 ), is a good measure to judge the performance of a given system as a qubit 2 ^. In 




(8) 



ih— = Hm 



(9) 
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actual physical systems, the leakage has a twofold origin, the driving used to switch between 
the basis states produces transitions to other levels besides the ones of interest, and the 
interaction with the environment. In this work we will only analyze the former without 
considering the possibility of ionization, so the electron remains bounded while the driving 
is applied. 

The material that follows is quite standard, nevertheless we include it for the sake of 
completeness. Writing the quantum state as a superposition of all the eigenstates 

^(f,t) = Y,c n (t)e- iE " t ^ n (r) (10) 

where 

H $ n (r) = E n <f> n (r), (11) 
is the eigenvalue problem, and replacing Equation [10] in the Schrodinger equation, [9J we get 

c n (t)e- lE ^ h H$ n = * ]T ^ {cn(t)e- iE ^ n ) S B . (12) 
which can be solved for Ck(t) using the orthogonality of the $ n 's^, 

.dc^}_ = S^ Cn{t){ $ klVext ^ n y. kn t_ (13) 

71=0 

Introducing the explicit for of the external driving, Equation [SJ we get that 

.dc k {t) 



A cos(u;t) Cn{t)ZkJ Wknt (14) 



at 

n=0 

where u kn = (E k - E n )/h, and 

Z kn = ($fc|z|$ n ), (15) 

clearly, the time-dependent probability that a given state, say k, is occupied at time t is 
given by \c k (t)\ 2 

Equation [TD shows that, even if the driving is in resonance with the frequency u res = 
(E q 2 — E q i)/h, there are transitions to all the bounded states allowed by the dipole selection 
rules so, unless somehow the matrix elements Z kn preclude this possibility, for large enough 
time all the terms in the superposition Equation [TU] will have non-negligible contributions. 
This fact, inevitably, produces a large and undesirable leakage. 

The points made above suggest that a proper design of the nano-structure would lead to 
negligible matrix elements Z kn for the unwanted transitions. This is the reason to choose 
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tips] 

FIG. 4. (Color on-line) Time dependent probability of finding the system in the lower state of the 
qubit (|cqi| 2 ) for different devices and for an rf field pulse of strength Aq at the resonant frequency 
uj res = (E q 2 — E q i)/h . The system is prepared in the lower qubit state for t = 0. The black solid line 
corresponds to the time dependent probability for the Device 1 and red dashed line corresponds to 
the time dependent probability for the Device 2. (a) Aq = 0.27 meV/nm, (b) Aq = 25.0 meV/nm 
and (c) Ao = 52.5 meV/nm. 

structures that single out a couple of eigenstates whose spatial extent is quite different 
from all the other ones since, clearly, this is a economical way to reduce the transitions 
to eigenstates that are not those of the qubit. The naive picture that says that the basis 
states of a qubit can be chosen as two well separated states of a physical system is really 
hard to be found, in particular when the requirements of fast enough operations also must 
be accomplished^. From a physical point of view, to obtain faster operation times it is 
necessary to draw on stronger external drivings. By two well separated states, it is meant 
that the energy difference between any other state with the qubit states is much larger than 
the energy difference between the qubit states. 

The matrix elements Z^ n can be obtained for all the bounded states of Devices 1 and 
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2 with very high precision so, the time evolution of the electronic quantum state results 
from the integration of Equation [TH The numerical integration was performed using high 
precision Runge-Kutta algorithms, taking into account all the bounded states that each 
device possess. 

Figure S] shows the time evolution of the occupation probability of one of the qubit basis 
states, the initial condition is ty(t = 0) = \q). The Figure shows the time evolution for three 
different driving strengths and two Devices, those whose eigenstates electronic densities are 
shown in Figure [3j Clearly, as the driving increases its value, the time evolution of the 
state becomes less and less harmonic. Anyway, for all the driving strengths shown, a fast 
switching between states can be easily achieved. The departure from a simple oscillatory 
behavior observed for larger driving strengths, and that |c q (t)| 2 < 1 for t > 0, shows that the 
driving is producing a superposition of many different eigenstates. On the other hand, from 
the three panels of Figure HJ the dependence of the switching time on the driving strength 
is manifest. 

As has been said above, that |c g (t)| 2 does not reach the unity for t > 0, is manifestation 
of the probability leakage, i.e. the state does not switch perfectly between the two qubit 
basis states. Anyway, for systems with a finite number of states, there is a finite probability 
that the state of the system returns to the initial state for a large enough evolution time. 
A potential well has indeed a finite number of eigenstates and, in many cases, the numer- 
ical integration of Equations [14] imposes a further reduction of the number of eigenstates 
effectively considered. For these reasons it is useful to introduce time-averaged quantities to 
qualify the dynamical behavior of the system. 

Together with the instantaneous leakage, 1 — (|c g i| 2 + c q2 \ 2 ), it is customary to introduce 
the time-averaged Leakage, L p , which is defined as 



where T is a large enough time that, in principle, can be taken equal to several periods of 
the external driving. Figure shows the behavior of the time-averaged Leakage for the two 
Devices previously defined as a function of the strength of the external driving. As can be 
appreciated, the time-averaged Leakage depends quadratically on the driving strength. On 
the other hand, Figure [5] shows that, despite that Device 2 differs from Device 1 in just one 
single eigenstate localized in the innermost potential well, the Leakage of both devices differ 




(16) 
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FIG. 5. (Color on-line) Time- averaged leakage L p for an rf field pulse of strength Aq at the resonant 
frequency oj res = (E q 2 — E q \ )/h as a function of the pulse strength for the Device 1 (black circle 
dots and line) and the Device 2 (red squared dots and line). 

in two orders of magnitude. It is worth to remark here that the spectrum of both devices 
are, essentially, the same (see Figure [2]). 

The good performance of Device 1 as a qubit can be further emphasized, looking at the 
dynamical behavior of the quantum state, when the system is forced with an off-resonance 
driving. 

Figure [6] shows the behavior of the Leakage as a function of the frequency of the external 
driving. It is assumed that the driving potential has the same properties that the one 
in Equation |HJ The Figure shows the data obtained for the two larger driving strengths 
showed in Figure H] since these two are the worse cases. As can be appreciated, for the 
smaller driving strength the driving frequency can be off-resonance up to ±10% without 
changing appreciably the Leakage. The size of this tolerance interval is closely related to 
the ratio between the eigen-energies differences with the driving strength. For the larger 
driving, the tolerance interval for frequencies smaller than ger as the 

tolerance interval for the smaller driving. The peaks in both sets of data signal that the 
transition probability to other states, which are not those of the qubit, gets bigger, spoiling 
the switching between the qubit states. 

It is worth to mention here that when the architecture of the quantum dot separates 
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FIG. 6. (Color on-line) The Leakage experimented by the Device 1 when the external driving is 
off-resonance. Note that both axis scale are normalized to the on-resonance values. The (red) 
triangle dots correspond to the a driving strength of Aq = 25.0 meV/nm and the (black) solid dots 
to Aq = 52.5 meV/nm, respectively. These driving strengths correspond to the cases b) and c) of 
Figure 0] 



two eigenstates, the residual leakage is produced by the small overlap between the qubit 
states and all the other eigenstates. Actually, because the small barrier that separates the 
potential wells, there is a non-negligible portion of the qubit states that lies in the outermost 
well. A better design would reduce or eliminate this portion enhancing the good behavior 
of the qubit. In the next Section we present a different model potential that possess exactly 
all the desirable properties that we have mentioned so far: separates in a potential well the 
two lowest lying eigenstates, and reduces the overlap between the qubit states and the other 
states to almost negligible values. Regrettably the model potential, and the parameters that 
characterize it, does not correspond to a nano-device already proposed. Anyway, since we 
have tested only a very limited set of options, it is possible that the discussion of the ideas 
presented here helps to find better designs for qubits based on multi-layered quantum dots. 
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FIG. 7. (Color on-line) a) The Leakage vs the driving strength Aq, the value of V is shown in the 
figure, while u> = 1 and 7 = 0.1. The black solid line corresponds to the device shown in panel 
b) and the dashed black line to the device shown in panel c). Note that the devices differ in the 
number of eigenstates localized in the innermost well, two and three states in panels b) and c), 
respectively. In panels b) and c) the black dashed line represents the decaying sinusoidal potential, 
in units of A. The color lines correspond to the electronic densities of several eigenstates. 



V. OTHER SYSTEMS 



So far, we have not been able to find an adequate set of parameters (radii, materials 
and so on) to design a quantum dot model that separates as qubit states eigenstates with 
n r = 0. The state with n r = i = has too many advantages to rule out too fast the search 
of a quantum dot such that one of the two separated eigenstates be it. By separated states 
we mean that two eigenstates are well localized in a given potential well, while all the other 
bounded eigenstates are localized in other potential wells. 

Anyway, the separation of the lowest lying eigenvalue, n r = £ = , and an excited state 
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with t = 1 can be achieved in the potential 



V c (r) = -V Q e~ 1T sin (wr), (17) 

where Vo, 7 and u are constants. The potential in Equation [T7] has some common properties 
with step-like potentials like the one in Equation The parameters of the potential and 
its shape are quite different from those commonly found in nano-devices, so in this Section 
we use atomic units. 

Figure [7J shows the Leakage suffered by the two "devices" depicted in panels b) and c) 
that has two and three states localized in the innermost well, respectively. The parameter 
that drives the change from two to three separated states is, in this case, V). The two lowest 
lying eigenvalues have quantum numbers (n r = 0, i = 0) and (n r = 0, 1 = 1). We call V C the 
critical value such that for Vo < V Q C there are only two well localized states in the innermost 
well. For Vo > V C there are more than two localized states in the innermost well, in our case 
V C « 2.09 

The eigenvalues and eigenstates of the one-electron Hamiltonian with potential [T7] were 
obtained approximately using the FEM-DVR method (finite-element method plus discrete 
variable representation), see^— . 

As can be seen in Figure [7J a) the Leakage of the device with three separated states is 
larger than the leakage suffered by the device with only two states separated. Again, the 
Leakage is a quadratic function of the driving strength, A . 

The influence of the number of separated states in the Leakage can be put further in 
evidence analyzing its behavior when Vo is swept from a smaller value than the critical to a 
larger value than V C . This behavior is shown in Figure [HJ As can be appreciated from panel 

a) there is a sudden change in the Leakage around the critical value, in fact the Leakage 
changes three orders of magnitude when Vq goes from smaller values than the critical to 
larger ones, irrespective of the driving strength, at least for the range of values analyzed. 
The jump can be rightly attributed to the presence of other states than those of the qubit 
since the others quantities involved (in the Leakage) are continuous, for example, Figure [8] 

b) shows the behavior of the resonance frequency as a function of Vo. 
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FIG. 8. (a) Time-averaged leakage L p for an rf field pulse of strength A = 0.001 (black solid dots 
and line) and Aq = 0.01 (triangle dots and dashed line), at the resonance frequency, as a function 
of Vq. (b) The resonance frequency as a function of Vq. 

VI. DISCUSSION AND CONCLUSIONS 

The ability of multi- wells and barriers spherical potentials to separate a subset of eigen- 
states in a potential well could be a tool to better designed nano-devices. This capacity 
has not been yet systematized, there are not exact results or theorems that allow a system- 
atic search of model potentials with all the desirable characteristics, despite it seems fairly 
general. 

As has been said above, once two eigenstates are separated in a potential well from all 
the other eigenstates, the residual leakage observed during the switching between them is 
attributable to the overlap between the qubit states and the other eigenstates. A higher 
barrier between the potential wells could reduce the overlap but, so far, we do not know of 
quantum dots built from three different semiconductor compounds. 

The advantage of the ground state as one of the qubit basis state is obvious, it can be 
pinpointed by cooling methods while other states require more sophisticated means to force 
the electron to actually occupy one of them. 

Our results contribute to show that there are, yet, a lot of improvements that can be 
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made to the design of qubits based on semiconductor quantum dots. The huge amount of 
materials and geometries offer ample possibilities to tackle the drawbacks that have marred 
the development of a reliable quantum dot based qubit. Of course the decoherence induced 
by the spin-orbit interaction, not considered in this work, stills remains as the heavier 
challenge. To minimize the effects of spin-orbit interaction the qubit states should be located 
in a potential well made of the material with the lowest spin-orbit interaction strength 
possible. 
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